Non-equilibrium transport theory of the singlet-triplet transition: perturbative 

approach 
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We use a simple iterative perturbation theory to study the singlet-triplet (ST) transition in lateral 
and vertical quantum dots, modeled by the non-equilibrium two-level Anderson model. To a great 
surprise, the region of stable perturbation theory extends to relatively strong interactions, and 
this simple approach is able to reproduce all experimentally-observed features of the ST transition, 
including the formation of a dip in the differential conductance of a lateral dot indicative of the 
two-stage Kondo effect, or the maximum in the linear conductance around the transition point. 
Choosing the right starting point to the perturbation theory is, however, crucial to obtain reliable 
and meaningful results. 
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I. INTRODUCTION 



With the evolution of mcsoscopic and nano-scale tech- 
nology and electron lithography, physicists gain more 
and more control over sub- micron structures and molecu- 
lar devices^— By exploiting cutting-edge semiconductor 
technology, experimentalists can now not only build ar- 
tificial atoms and molecules to trap and manipulate indi- 
vidual electrons and their spin^i 6 - but they are also able 
contact and gate real molecules; 7 - and build simple molec- 
ular devices out of them such as the single electron tran- 
sistor. These devices are promising candidates for our 
future electronics, and might possibly be used as build- 
ing blocks of spintronics devices and spin-based quantum 
computers, too£ Moreover, the unprecedented control 
of these minuscule structures opens new and fascinat- 
ing possibilities to build and study hybrid structures^— 
entangle electron spins^ 4 - and also to realize and study 
simple quantum systems in the close vicinity of quantum 
phase transitions under non-equlibrium conditions.— ~— 

Understanding the physical properties of these tiny 
electronic circuits represents a major challenge to theo- 
retical as well as to experimental physicists. Due to their 
extremely small size, electron-electron interaction is of- 
ten dominant in these structures, and moreover, as men- 
tioned before, they are operated under non-equilibrium 
conditions. Although theoretical physicists devoted a lot 
of effort to design methods that are able to tackle these 
difficulties, so far none of the proposed methods proved 
to be entirely successful in describing strongly inter- 
acting non-equilibrium systems: Monte Carlo methods 
are presently unable to reach the required precision ; 20 ! 21 
Bethe Ansatz methods can be used for a few models only, 
and are still in an experimental stagej22r— and pertur- 
bative renormalization group methods can only reach a 
particular region of the parameter space.— ~— Maybe nu- 
merical renormalization group methods arc currently the 
most reliable techniques to study these non-equilibrium 



systemS ) 24 i 29 ' 30 however, they scale very badly with the 
number of states involved, and to compute the trans- 
port through just two levels in the presence of interaction 
seems to be numerically too demanding. 

In view of the above-described situation, pertur- 
bative methods such as iterative perturbation theory 
(IPT) (2Ir— e.g., are of great value, especially, since many 
well-contacted or almost open systems are de facto in 
the perturbative regime. For this reason, several molec- 
ular electronics groups combine the GW method with 
standard ab initio tools to compute transport through 
molecules^ Unfortunately, however, very little is known 
about the reliability of this approach. It is, e.g., rather 
easy to obtain spontaneous local symmetry breaking 
(e.g., magnetic moment formation), and thereby obtain 
physically incorrect results^ 

In the present paper we investigate systematically the 
perturbative approach for the simplest possible non- 
trivial molecular system, a two-level Anderson model 
with two degenerate or nearly degenerate electron levels. 
This simple model describes many experimental systems, 
and it displays rich physical properties. Having an odd 
number of electrons on the two levels, e.g., one can re- 
cover an almost SU(4) symmetrical Kondo stated - — In 
case of two electrons, on the other hand, one finds a tran- 
sition (cross-over) between a state, where the two elec- 
trons are bound into a singlet, and another state, where 
the two electrons' spin is aligned into a triplet, which is 
then completely or partially screened by the spin of the 
conduction electrons in the leads] 16 ' 17 ' 41 " — Here we fo- 
cus our attention to this second, most interesting regime, 
also nicknamed as the singlet-triplet (ST) transition. We 
focus on the most realistic case^ where two indepen- 
dent conduction electron channels are coupled to the dot 
levels, and therefore the singlet-triplet transition is just 
a cross-over 4Ir— This must be contrasted to the case of 
a single coupled conduction electron channel, when the 
transition is a true Kosterlitz-Thouless quantum phase 
transition ! 16 ' 17 Non-equilibrium transport in this latter 
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case has been investigated recently using the so-called 
non-crossing approximation (NCA)4^ 

Describing the case of two coupled channels under 
non-equilibrium conditions represents a theoretical chal- 
lange. The deep singlet side of the transition is acces- 
sible through perturbative renormalization group meth- 
ods, which were used successfully to analyze the non- 
equilibrium transport in this perturbative rcgime4£ Here 
we focus on the transition, and investigate, whether a 
simple perturbative approach is able to account for the 
transition between the two correlated states. To our sur- 
prise, perturbation theory captures all qualitative fea- 
tures of the transition, including the two-stage Kondo ef- 
fect, or the conductance maximum around the transition 
point4^~— However, to obtain accurate perturbative re- 
sults, one needs to expand around the best possible non- 
interacting theory.— ~— This we achieve by introducing 
self-consistently determined countertems, which is essen- 
tially equivalent to doing perturbation theory around the 
Hartrcc solution. 

The simple perturbative calculation is, however, un- 
able to capture correctly the Kondo scales, i.e. the en- 
ergy scales below which the electron spins on the dot 
are screened or bound to a singlet. As we show in a 
Ref. this shortcoming of the iterative perturbation 
theory can be overcome to a large extent by the so-called 
Fluctuation Exchange Approximation (FLEX) extended 
to non-equilibrium situations. 

The paper is organized as follows. In Sees. Ill Al and 
IIIB1 we introduce the non-equilibrium two-level Ander- 
son model, and describe the iterative perturbation the- 
ory used to solve the non-equilibrium Anderson model. 
In Sec. Ill CI we show how to obtain the transport prop- 
erties. Some computational details and the limitations 
of the iterative perturbation scheme are given in sec- 
tions III Dl and III El In Sec. IIII Al we discuss the results 
obtained for completely symmetrical quantum dots with 
equal level widths, while in Sec. IIII Bl rcsults for dots with 
more generic parameters are presented. Our conclusions 
are summarized conclude in Sec. IIVI and some technical 
details are given in Appendix lAl 

II. THEORETICAL FRAMEWORK 
A. Model 

The physical properties of a quantum dot with even 
number of electrons are usually well-captured by a simple 
two-level Anderson Hamiltonian. This can be written 
as a sum of a non-interacting and an interacting term: 
H = Hq + -Hint- The non- interacting part, Hq, can be 
further divided onto three terms, 

H = -ffcond + #hyb + -ffo.dot- (1) 

Here i? C ond describes conduction electrons in the leads, 
-ffhyb accounts for the hybridization between the con- 
duction electrons and the electrons residing on the dot, 



and -ffo,dot stands for the non-interacting part of the dot 
Hamiltonian, 

#0,dot = ^Eid^dur , (2) 

I .o 

with <f i(J the creation operator of a dot electron of spin a 
on dot level i of energy £j (i e (+, — )). 

The precise form of H con d and -f/hyb depends on the 
geometrical arrangement of the quantum dot. In case of 
a lateral quantum dot close to pinch-off, there is only a 
single conduction electron mode coupled to the dot, and 
correspondingly, 

^cond = E ^ aC \aa C £,a° ' 

^yb = E t<*i{c\ aa di° + h.c). (4) 

Here c^ aa creates a conduction electron in the left or 
right lead (a 6 (L,R)) with energy £ Q = £ + fj, a , with 
fi a = eV a the bias applied on lead a. The t a i's denote 
tunneling matrix elements between lead a and dot level 
i, and for the sake of simplicity here we assume them to 
be spin- and energy independent. 

Vertical quantum dots^r— or nanotubes^i are, on the 
other hand, often connected to leads with a large surface 
and therefore in each lead there are typically many differ- 
ent conduction electron modes which couple to a given 
dot state. However, for each dot state i and lead a, 
one can construct a single linear combination of modes, 
d\ a — > c| iaa , which hybridizes with di a . In many cases, 
one can assume (based on symmetry considerations or in 
view of the chaotic shape of the dot wave functions and 
the large number of lead modes) that i aa arc indepen- 
dent of each other for i = ± i.e. {c| + a a , cj ; _ ;Q , jCr } = 0. 
Under this assumption, H con d and -f/hyb assume the form, 

-^cond = £,a C \ia<r C i ia ° ' 

H lyh = E t ai (cl ia A<r + h.c.) • (6) 

Notice that the non-equilibrium condition is taken into 
account through a lead-dependent shift of the chemical 
potentials, and the interaction between the conduction 
electrons is neglected. 

In general, the hybridization matrix elements t a i are 
independent of each-other. However, for simplicity, in 
this paper we focus on completely symmetrical quantum 
dots and assume further that one of the dot levels is even 
(+), while the other is odd under reflection (— ). As a con- 
sequence, for a lateral dot the tunneling matrix elements 
have a simple structure: th,+ = tn,+ and i^,- = —tn,-. 
Although this choice may seem to be too special, this 
symmetrical model still captures most of the generic fea- 
tures of a two-level quantum dot. 43 ' 44 
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FIG. 1: The antisymmetrized vertex function containing the 
Coulomb interaction and the Hund's rule coupling. 



In our model, electron-electron interaction on the dot 
is taken into account by the term: 



U 



J s 2 , 



(7) 



relations, 



U 3J N 
2 + T 



J2 d L d L d i<t d ii 



U J 

1 ~ 4 

U J 

2+4 



Yj 4a d ]a d jS d t. 



(10) 



where a = — a. To simplify perturbation theory, it is 
useful to collect all fourth order terms and rewrite : H m t ■ 
in terms of the antisymmetrized interaction vertex shown 
in Fig.IJr^r'.as 



Hint J ^{a dl„,dl„~r 



no 3<r ma' 



(n) 



z,j,m,n 
ct,(t' a .a 



where 



dl^dia and S 



\ J2 d i<r°(T<?'dia< With 



a denoting the Pauli matrices, a = (a x ,a y ,a z ). Here 
the coupling constant U denotes the on-site Coulomb in- 
teraction, and accounts for the charging energy of the 
dot, while J stands for the Hund's rule coupling, favor- 
ing a ferromagnetic alignment of the dot electron spins 
(J > 0). In this work, we focus our attention to the 
regime where approximately two electrons reside on the 
dot. This assumption is implicit in the way we expressed 
the Coulomb term in Eq. (J7J, which obviously possesses 
a particle-hole symmetry discussed later. 

Although Eq. (J7]) provides a rather symmetrical and 
compact expression for the interaction term, for the pur- 
pose of a systematic perturbation theory, it is useful to 
rewrite it in a normal ordered form, 



Using the fact that r^ n ™ CT is anti-symmetrical under 
exchanging ia o na and ja' -O- ma' , it is easy to read 
out the non- vanishing matrix elements of T?" n ™ CT from 



Eq. (|T0)) . Although the vertex T has formally 256 ele- 
ments, as a consequence of spin and level-conservation in 
Eq. (|10j) . only thirty- two of these enter actual calcula- 
tions. The advantage of the symmetrical form, Eq. 
is that it allows one to treat the Coulomb interaction 
and the Hund's rule coupling on equal footing: one can 
formally (and numerically) evaluate the self energy di- 
agrams shown schematically in Fig. [21 with the specific 
form of the interaction hidden in the internal structure 
of these diagrams showing in Fig. Q] 

Before proceeding with the presentation of our non- 
equilibrium perturbation theory, let us shortly discuss 
the symmetry properties of the Hamiltonian. Clearly, 



Hi n t ='■ H, nt : — I H — — 



W_ 
2~ 



3J 
4 



i . a 



(8) 



and treat only the normal ordered part as a perturbation, 
: H m t while adding the second term in this expression 

to -ffo.dot: 



a) 



b) 




3J 3U 
a -> Ei = ei-— — 



(9) 



We remark that for a two-level dot close to half filling, 
6j ~ 0, and the point e+ = e_ = corresponds to the case 
where the two levels are exactly degenerate, and there are 
exactly two electrons on the dot ((^2 ia Tii a ) = 2). 

For later purposes, let us give here the explicit form of 
: Hi n t '■ as obtained by exploiting fcrmionic commutation 




FIG. 2: Countcrterm (a), first order (b) and second order 
(c) contributions to the self-energy. The lines denote non- 
interacting Keldysh Green's functions. 
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H is invariant under SU(2) spin rotations. This implies 
that, in the absence of an external magnetic field, all 
Green's functions are diagonal in the spin labels and also 
independent of them. 

In addition to spin SU(2) symmetry, depending on the 
specific Hamiltonian parameters, H may also possess a 
discrete electron-hole symmetry. In fact, throughout the 
present paper we shall restrict ourself to the simple case, 
where the two levels e± are shifted symmetrically away 
from 0, 



e± = ±A/2 



(12) 



Then, for |i a +| = |t a — |j the Hamiltonian is also invariant 
under the following transformation: 



(4 ifl d\ ti ) o (d_ ti , - d_ it ) 



l) (- C " 



(13) 



This symmetry implies that the occupation of the dot is 
exactly (Y). 7i«r) = 2 for any interaction strength and 
splitting A, and also implies the relation 

between the spectral functions g± of the two dot levels. 
[For a precise definition, see Eq. ([33])] 

Finally, at the special point, e + = e_ = 0, wc have yet 
another electron- hole symmetry, satisfied even for \t a+ 1 ^ 



(<4,t' 4,4.) ( d «,U - d «,t) , 



(14) 



This second symmetry impleas the spectral function re- 
lations, g+(u>) = g + (—u) and Q—(oj) = Q-(—u>). 



B. Out of equilibrium iterative perturbation theory 

To describe the transport properties of a quantum dot 
under out of equilibrium conditions, we shall apply the 
so-called Keldysh Green function technique^ However, 
before doing so, let us further reorganize our Hamilto- 
nian. Clearly, interactions can substantially shift the ef- 
fective values of the dot energies. Therefore, to account 
for this trivial but possibly large renormalization effect, 
we apply a counterterm procedure: We take as a non- 
interacting part the following Hamiltonian 



Ha — H. 



cond 



-ffhyb + ^ £ia d\ a d\ c 



(15) 



while the rest of the Hamiltonian is treated as a pertur- 
bation, 

-Hint = : Hixft : +-ffcount ) (16) 

tfcount = ^2(Ei - s ia )dl a d ilT . (17) 



In other words, we expand around fictitious (effective) 
dot levels, e^, at the price of treating the counterterm 
-ffcount also as a perturbation. This approach is essen- 
tially the same as performing perturbation theory around 
the Hartree theory, and allows us to extend the range of 
validity of our perturbation theory by making a proper 
choice oi Ei a (see below). 

Our primary goal is to compute the non-equilibrium 
Keldysh Green's functions, G K ' K , with the index pair 
(ft, ft') corresponding to the branches of the Keldysh con- 
tour. With this notation G 11 stands for time-ordered 
(G T ), G 1 ' 2 for lesser (G<), G 21 for larger (G>), and G 2 ' 2 
for the anti-time-ordered (G T ) Green's functions. These 
are related to the usual retarded (R), advanced (A) and 
Keldysh (K) Green's functions by the usual relations^ 2 - 



G R = G T - G< 
G A = G T - G> 

G K = G T + G f 



(18) 
(19) 

(20) 



and also satisfy the restriction, G T + G T = G < + G > . 

To determine the interacting dot Green's functions, 
G1^ Va i , we perform perturbation theory in H int by tak- 
ing Ho as a non-interacting Hamiltonian. The Green's 
function can thus be expressed in terms of the self-energy 
VrT , and the unperturbed Green's functions, g^. Va , by 
the Dyson equation, 



G(w)- 1 = g(w)- 1 - E(w) 



(21) 



where a matrix notation A^ -,, — > A has been in- 
troduced. The non-interacting Green's functions corre- 
sponding to the resonant level model Hq can be easily 
computed, and are listed in Appendix lAl Up to second 
order in the interaction, the self-energy can formally be 
written as 



S(w) = £( count > + £W + S (2)( a 



(22) 



where the different terms denote the counterterm, the 
first and second order terms shown in Fig. [21 An explicit 
evaluation of the diagrams in Fig. [2] gives the following 
results: 



vj(i) K "' =s , s 



\ pjc rn<7 pi a n <j 
SkSs / j 1 ia tut'" L j'a' ma'- 

j' ,rn' ,n' 



duii f du>2 



2?T J 2tt 

— oo 



x g^ a , (wi) g^ m , a „ (u> 2 ) K - W2 - w) , (23) 

where the Keldysh sign s K has been introduced to ac- 
count for the sign change of the interaction term on the 
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Kcldysh contour: si = 1, S2 = — 1. In these expressions 
we assumed that the z-component of the spin is con- 
served, and correspondingly, the Green's functions and 
the self-energies are diagonal in a. 

In principle, Eqs. ([2"T]). (|22|) . and <(23j) give a com- 
plete perturbative description of the quantum dot. How- 
ever, they depend paramctrically on the so far unspeci- 
fied levels, e± a - Following the concept of iterative per- 
turbation theory, we determine these quantities selfcon- 
sistently from Eqs. (|22p and (|2"Tj) by requiring that g and 
G give the same occupation numbers 
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(24) 



Condition (|24|) can be fulfilled by tuning the effective 
levels. In spite of the simplicity of this "iterative pertur- 
bation theory" , the results obtained with it show remark- 
able agreement with the ones obtained with more ad- 
vanced techniques for several equilibrium systems.— Un- 
fortunately, a well established and exact non-equilibrium 
impurity solver being not available yet, one can only 
judge the validity of this perturbative approach by pos- 
sibly comparing with results obtained by some reliable 
equilibrium methods, and by investigating the internal 
consistency of its. As discussed in Section III El the range 
of applicability turns out to be limited to small to mod- 
erate values of U and J. 

In principle, Eq. (f24| could and should be solved for fi- 
nite biases. However, under non-equilibrium conditions, 
the off-diagonal elements (d\ a d,2cr) become also finite, 
and additional constraints related to current conservation 
may also emerge. In this work, we rather use a different 
and simple strategy: we determine ii in equilibrium, and 
then fix it also for the non-equilibrium calculations. For 
not very large bias voltages this leads to a stable solution 
of the self-consistency equations. 



C. Transport properties 

Having the Green's functions G at hand, we can calcu- 
late transport properties such as current and conductance 
based upon the Meir-Wingreen formula^, 



2h^ 



ckuTr 



(r i -r«)G<H + 



— 00 

+ (/ L ( w )r L -.feHr«)(G«H-G^H)j .(25) 

Here f a (oj) = f(u) — fi a ) is the Fermi function in lead a, 
and the Green's functions arc matrices in the level indices 



only. The difference between the formulas for lateral and 
vertical dots appears only in the different structure of the 
T-matriccs defined as 



(IXOlat = 27T N t ai t* 
in the lateral case and 



(26) 



(Iy)vert =Sij 27T N t ai t* aj , (27) 

for a vertical quantum dot. These are related to the 
width r± of levels i = ± through 



In our case, the ground state of the system is a Fermi 
liquid^. Therefore, according to Nozieres's Fermi liquid 
theory, at T = temperature (and V = bias), quasi- 
particles at the Fermi energy scatter elastically from the 
impurity, and their scattering process can be described in 
terms of simple phase shifts. Since our model is invariant 
under reflection for V = 0, and since one of the levels 
is assumed to be even while the other one odd, we can 
characterize the scattering process at the Fermi energy 
by just two phase shifts, Si, associated with the two dot 
levels. These phase shifts are related to the retarded 
equilibrium Green's function through the Fermi liquid 
relation^ 



7T 

2 



arctan 



fImGt> 



0,F = 0) 



UeG£> = 0,V = 0) 



(28) 



In the limit of infinite bandwidth, the phase shifts are 
also related to the occupation of the levels di a by the 
Friedel sum rule^ 



2- - (m) 

7T 



(29) 



where (rii) = X)o-( n «cr)- However, for a finite conduction 
band cut-off Eq. (|29[) is only approximate, and Eq. (|28|) 
gives a more reliable way to determine the phase shift. 

The T = linear conductance is directly related to 
the phase shifts above through the Landauer-Biittiker 
formula. For a lateral dot, an elementary calculation 

A 1 

give: 



.43 



G 



lal 



2e 2 



sin 2 (£4 



5-] 



while for a vertical dot one obtains^ 

2e 2 



/Tver 



h 



(sin 8+ + sin S- , 



(30) 



(31) 



These equations are very instructive, and help us to un- 
derstand the transport properties of the dots. They im- 
ply, e.g., that at complete resonance, 6+ — 5- = n/2, 
the conductance of a vertical dot is 4e 2 //i, while that of 
a lateral dot vanishes due to the interference of scattering 
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states. These interference effects are expected to show up 
also in the non-equilibrium results. 

The equations above are also useful to test our nu- 
merical calculations. The linear conductance values in 
Eq. (|30|) and in Eq. ([31]) can be compared to the equi- 
librium limits of the differential conductances obtained 
by a numerical differentiation of the current as a func- 
tion of the bias voltage dI/dV\v=0: as computed from 
the Meir-Wingreen formula. 

D. Computational details 

In the numerical calculations wc represented the 
Green's functions using a finite uniform mesh of N fre- 
quency points in the range —Q/2 < lo < fl/2. Depending 
on the accuracy needed, we used slightly different values 
for O and N. To study the range of applicability of our 
method, we used N = 2 15 points and il = 500C7. To 
compute the spectral functions, we used N = 2 17 and 
Q = 1000J7, while for the current and differential con- 
ductance plots we chose N = 2 17 and ft = 200U. All 
results presented in this paper were obtained at T = 
temperature. 

The advantage of the iterative perturbation theory is 
that the calculation of the second order self-energy is only 
seemingly cumbersome. Rewriting Eq. (|23l) in time do- 
main, one obtains a simple multiplication of the various 
components of the Green's functions, 

y(2) KK ' /.\ _ \p pjV mu" r iVnV" 

^ ii'<j\ b ) ~ 5 kV L ia na'" L j' a' ma" A 

o-VV" 

x S$* (*) <&tv" (*) fl£w« (-*) • (32) 

Thus the self-energy can be computed by just perform- 
ing a Fast Fourier Transformation to get g(t), evaluating 
E(t) in the time domain, and then transforming it back 
to frequency space. 

As already mentioned before, for any given set of pa- 
rameters we first carried out a calculation in the absence 
of voltage, V — 0, and determined the levels, e±. Then 
we computed S(w) and G(t) using Eq. (|52")| , computed 
the spectral functions 

Q ia (u,) = --lmG% a (uj), (33) 

7T 

as well as the total spectral function, Qt{u) = 
J2ia fto-(w), and finally evaluated the current using the 
Meir-Wingreen formula, Eq. (|2"5|) . The differential con- 
ductance has been obtained by direct numerical differen- 
tiation of the current, G(V) = dI(V)/dV. 

E. Limitations 

Perturbation theory has a limited range of validity. For 
large values of the interaction, U, a spontaneous symme- 
try breaking occurs, whereby the occupation of the dot 
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FIG. 3: (Color online) Hysteresis of the occupation number 
as a function of level splitting (A) for J/U = 0.0 and T±/U = 
T/U = 0.4. 

levels i = ± becomes unequal even for A = e + — e_ = 0, 
(n + ) =/= (ri-). This phenomenon is similar to the sponta- 
neous moment formation found by Anderson , 35 ' 55 which 
wc eliminated here by setting the Green's functions ex- 
plicitly proportional to ~ S aa i . In a perturbative calcu- 
lation one should always stay away from these regions in 
parameter space to avoid spurious solutions. Fortunately, 
as shown in Fig. [3] we can easily identify these " danger- 
ous" regions by computing the occupation numbers (n±) 
as a function of level splitting, A/U, since {n±) exhibit 
hysteresis there. 

Carrying out a similar hysteresis analysis for each value 
of the parameters U/T and J/T, we obtain the stability 
diagram, Fig. [H delineating the region of applicability of 
our perturbative approach (for simplicity, here we took 
r± = r). Clearly, one needs to keep both U/T and J/T 
moderate to stay within the range of applicability of it- 
erative perturbation theory. 



III. RESULTS AND DISCUSSION 

A. The symmetric case (r+ = T_) 

The parameter space of the two-level quantum dot is 
huge. Therefore, to get more intuition, we first restrict 
somewhat the number of independent parameters by as- 
suming that both levels couple to the leads with equal 
strength, r± = T. We assume further that the dot is at 
the particle-hole symmetrical point, e+ = — e_ = A/2, 
where the number of electrons on the dot is fixed to 
n = CT ( n io-) = 2. Before presenting the numerical 
data, let us briefly discuss the structure of the ground 
state as a function of J, U, and A in case of small tun- 
nclings, T. 

In the absence of tunneling to the leads, T = 0, the 
states of the dot can be grouped into three singlets 
(with spin S = 0) and a triplet (S = 1). Two singlet 
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FIG. 4: Stability region of perturbation theory for V = F±. 
Empty circles mark the region of hysteresis, while filled circles 
indicate the stable region, where PT gives meaningful results. 
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FIG. 5: Evolution of the singlet and triplet states of the 
isolated dot for fixed J > 0, as a function of A. Energies 
are measured from the energy E t of the triplet state. The 
singlet-triplet transition point is also indicated (S-T). 



states, \s±) correspond to putting two electrons on the 

upper/lower level, while a third state, |s_| }, is a singlet 

with one electron residing on each level. The rest of the 
two-electron states are spanned by a triplet state, \t, m) 
(in = 0, ±), where the spin of the electron on level i = + 
is aligned with that on level i = — . Fig. [5] shows the evo- 
lution of these four levels as a function of the splitting 
A for finite J. For A = the triplet state is the ground 
state, and the three singlet states form degenerate ex- 
cited states of energy 2J. Increasing A, however, pulls 
down the state |s_), which becomes the ground state for 
A > 2J. 

This picture is slightly modified if the dot is coupled 
to the leads. In this case, for generic couplingS ) 43 i 56 
the transition becomes a smooth cross-over, so that the 



triplet and singlet states are adiabatically connected. 
However, the various regimes are described by rather dif- 
ferent physical pictures. For large values of A, the two 
dot electrons form a local singlet on the dot, while for 
J <C |A| they are aligned into a triplet, which then cou- 
ples to the leads through an exchange coupling, and is 
screened by a two-stage Kondo effect for T + < T_ at 
some temperatures < T^. In this regard, the case 
r_ = T + is rather special, since then only a single Kondo 
scale appears, Tj£ = As one approaches the tran- 
sition region, A rj 2J, from the triplet side, the Kondo 
temperatures T K increase, until the singlet state s_ and 
the triplet states are all inseparably mixed into a sin- 
gle quantum state, characterized by a Fermi liquid scale, 
T££L££ These three regions are sketched in Fig. O 



1. Spectral functions 

To gain insight to the effects governing the conduc- 
tance through the dot, it is useful to study first the equi- 
librium spectral functions. 

Case of no Hund's rule coupling, J = 0. Let us first 
investigate in the absence of Hund's rule coupling, J = 0. 
In this case, there is no triplet region, rather, changing 
A, one finds a transition between the two singlet states 
|s_) — > |s_|_) through a peculiar Kondo state formed by 
all three singlet states and the triplet (see Fig. [S]). Fig. [7] 
shows the total spectral function, qt(u>), in this case, for 
T/U = 0.565, as a function of level splitting, A/U. 

For A/U ~ and J/U = 0, the ground state of an 
isolated dot is sixfold degenerate. When the quantum 
dot is connected to the leads, a strong Kondo resonance 
driven by the quantum fluctuations of this sixfold degen- 
erate state appears at cj = 0. The width of this peak 
can be identified as the Kondo temperature, T£. In ad- 
dition, one observes two shoulders at cj « ±U, which can 
be identified as the Hubbard peaks. 

Increasing A first slightly suppresses the Kondo res- 
onance at u) = 0, and finally, for A > splits the 
resonance into two side-peaks. These side-peaks can 
be identified as the singlet-triplet excitation peaks, and 
their position can be precisely determined from the level- 
projected spectral functions, Q±(uj). As shown in the 
inset of Fig. for large values of A they are indeed 
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FIG. 6: Sketch of the singlet, triplet, and transition regimes, 
as a function of A/U for a fixed J > 0. 
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FIG. 7: Total equilibrium spectral functions, qt(w), for 
J/17 = and T±/U = V/U = 0.565, and different values 
of the level splitting, A/U . Inset: Evolution of the peak posi- 
tions in q±(ui) as a function of A/U. The dashed lines indicate 
uj = ±A 



FIG. 9: Linear conductance of lateral and vertical dots for 
J/U = 0.2, from numerical derivation (ND) of the current 
at zero bias, and from the Fermi liquid relations (FL). In- 
set: Comparison of occupation numbers and phase shifts, S/tt. 
The conductance of a symmetrical lateral dot shows a maxi- 
mum of height 2e 2 /h around the ST transition. 



located at w « ±A, although the splitting takes place 
rather abruptly at A w 0.4, and the A— dependence of 
the maxima of g±(u>) is very non-linear. 

Finite Hund's rule coupling, J > 0. The effect of 
the level splitting A on the equilibrium spectral func- 
tion is presented in Fig. [5] for parameters J/U = 0.2 
and T/U = 0.565. The spectral functions behave quite 
similarly to those calculated for J = 0: A Kondo reso- 
nance appears for A = in qt{oS), which is gradually 
suppressed as A increases, and finally splits into two res- 
onances for large values of A. However, in this case, we 
should interpret the Kondo resonance slightly differently 
than before. For A = the ground state of the iso- 



lated dot would be a triplet state, and therefore turning 
on a small T, this spin triplet gets screened by the con- 
duction electrons in the leads through a triplet Kondo 
effect— Although T is relatively large in our case and 
comparable to the singlet triplet splitting, T ~ 2 J, nev- 
ertheless, this interpretation still holds, since the width of 
the central resonance (Tr-) is still smaller than the split- 
ting AE st = | A — 2J|. In fact, the resonance splits ap- 
proximately where AE st ~ Tk, and the positions of the 
split resonances follow approximately the straight lines 
uj w ±(A - 2J) for large A's (see the inset of Fig. |8j), in 
complete agreement with the ST transition picture. 




FIG. 8: Total equilibrium spectral functions, qt(u), for 
J/U = 0.2, V±/U = V/U = 0.565, for different values of level 
splitting, A/U. Inset: Peak positions of g±(ui) as a function 
of A/U. The dashed lines indicate cj = ±(A — 2J). 



2. Conductance 

Having investigated the equilibrium spectral functions, 
let us now turn to the discussion of non-equilibrium 
transport through the dot. We first focus on the linear 
conductance, Gy ln = dl /dV\v=o- This is plotted in Fig. [9] 
as a function of A both for lateral and vertical dots with 
a finite Hund's rule coupling, J/U = 0.2. We computed 
Gii n , on one hand, by numerically differentiating the cur- 
rent, Eq. (|25p. but we also computed it by extracting 
the phase shifts from the retarded Green's functions and 
the occupation numbers, and then using the Landaucr- 
Biittiker formula (Eq. ([30j for lateral and Eq. (JHU> for 
vertical dots). As shown in Fig. [9l the two procedures 
give identical results within numerical accuracy. 

Vertical and lateral dots exhibit very different charac- 
teristics. The linear conductance of a lateral dot shows 
a maximum around the ST transition and it is small 
on both sides of the transition, while the conductance 
of a vertical dot crosses over smoothly from a conduc- 
tance of Gii n w Ae 2 /h at A = to a small value for 
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FIG. 10: Differential conductance, G(V), for J/U = 0.2 and 
T±/U = V/U = 0.565 for different level splittings, A/U both 
for lateral (a) and for vertical (b) quantum dots. For vertical 
dots the suppression and the splitting of the Kondo resonance 
are clearly visible in G(V). For lateral dots, interference ef- 
fects mask the signal on the triplet side. 



A 3> 2 J. This behavior can be understood from Eqs. ([50)1 
and (f3lj). and the behavior of the phase shifts, shown 
in the inset of Fig. |9j In the triplet regime, A w 0, 
two conduction electrons are needed to screen the local 
spin S = 1. Therefore, both phase shifts are close to 
S± « 7r/2 by the Friedel sum ruleJ^ While for a ver- 
tical dot the contributions of the channels i = ± add 
up to the conductance, and amount in a total conduc- 
tance of Gun = 4e 2 //i, for lateral dots there is a destruc- 
tive interference (see Eg. I50|) . which finally amounts in a 
complete back-reflection of electrons and no conductance, 
G ~ i 44 i 49 i 56 Increasing A, the phase shifts gradually 
cross over to values, 5- ~ n and S + ps 0. This gives rise 
to a maximal conductance of G\ ln = 2e 2 /h somewhere in 
the vicinity of the ST transition for a lateral dot^ while 
it results in the monotonous decrease of the conductance 
of a vertical dot^ as shown in Fig. [9] The Friedel sum 
rule, Eq. ([2l?|) . is verified in the inset of Fig. |H1 where we 
compare the phase shifts to the occupation numbers. 

The differential conductance, G(V), can be observed in 
Fig. fTDJ For a vertical dot, the behavior of G(V) follows 
the naive expectations based upon the spectral functions' 
structure and the excitation spectra of the dots: For 
A = a clear Kondo resonance is seen at V w 0, and the 
singlet excitations give a broad resonance at V ps 2J/e. 
Increasing A, the central resonance is suppressed in am- 



plitude (as predicted by the Friedel sum rule) , but it also 
broadens somewhat as one approaches the ST transition 
point, where quantum fluctuations give rise to a some- 
what higher Kondo temperature, T^- i i 58 The excitation 
peak at V ps (2 J — A)/e shifts to lower energies and 
gradually merges with the central resonance. Finally, for 
A > 2 J the quantum dot is in the singlet region. There 
the central resonance splits up and a resonance is ob- 
served at V w (A - 2J)/e. 

While it is easy to understand the differential conduc- 
tance of a vertical dot, the conductance of a lateral dot is 
counter-intuitive. For A = there is a complete destruc- 
tive interference, and G(V) = 0. We emphasize that this 
is only valid for the symmetrical case, T + = F_, studied 
in this section. Increasing A, G(V) becomes non-zero, 
and the conductance gradually increases. The shape of 
G(V), however, does not have a simple intuitive expla- 
nation on the triplet side of the transition, as a result 
of the cancellation of the various contributions. Around 
the transition point and beyond that (A > 2J), however, 
the phase shifts are far away from 7r/2, interference ef- 
fects are suppressed, and the conductance shows features 
which arc quite similar to those found for vertical dots, 
discussed above. 



B. The asymmetric case, 7^ T_ 

In the previous subsection, we analyzed the case of 
equal tunneling rates, T + = T_. While this assump- 
tion may be a good approximation for some systems (for 
carbon nanotubes, e.g.);* it is violated for most quan- 
tum dots. 41 ' 42 Let us therefore investigate in this sub- 
section the generic case of unequal couplings, T + ^ r_. 
The range of applicability of our perturbative method is 
shown in Fig. [TTJ For simplicity, we fixed T + = T_/2 
throughout this subsection, and characterized the hy- 
bridization strength in terms of the geometric mean of 
r±, r = v/TVTT. A comparison of Fig. QI] and Fig. H 
shows that the range of applicability of perturbation the- 
ory is only weakly modified by the assymetry in the tun- 
neling rates. 

For r + 7^ r_, a few important differences ap- 
pear though compared to the case of equal tunneling 
strengths. When T + 7^ r_, the particle-hole symmetry, 
Eq. (p~3|) . of the Hamiltonian is not valid anymore. As a 
consequence, the number of electrons on the dot is not 
exactly two and the spectral functions also violate the 
electron- hole symmetry relation, and £>+(w) Q-(-uj). 
Electron-hole symmetry is only maintained for A = 0, 
where Q+(u>) = q + (—uj) and Q—(u>) = Q-{—uj) follow 
from the second electron-hole symmetry transformation, 
Eq. (fH]l . The second major difference is that, due to 
the two different coupling strengths, two separate Kondo 
scales emerge for the even and odd channels on the triplet 
side, T±4i 

Both differences are clearly visible in the spectral func- 
tions, shown in Fig. 1121 where the level-projected spectral 
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FIG. 11: Range of applicability of perturbation theory for 
r+ T_ , as extracted from the hysteresis of occupation num- 
bers. Filled circles indicate the region of stable perturbation 
theory. Couplings are measured in units of Y = vT+r_ . 



functions are also displayed. Clearly, only the A = 
spectral functions are electron-hole symmetrical, and 
even there, the spectral functions g + and Q- have a cen- 
tral Kondo resonance of different width. 

Apart from these differences, the overall evolution of 
the spectral functions as well as that of the phase shifts 
(see inset of Fig. [T3|) is quite similar to the one observed 
for r + = r_. As a result, the behavior of the linear 
conductance as a function of the level splitting, A, is 
also quite similar to the one obtained for T + = T_ (see 
Fig. EJl. 

The differential conductance is, however, quite differ- 
ent, and shows a much richer structure on the triplet side 
in this case. For A = 0, both phase shifts are pinned to 
5± = 7r/2 by electron- hole symmetry. The physical rea- 
son for this is, of course the formation of a Fermi liquid 
state through a Kondo effect in both channels on the 
triplet side. As a consequence, by Eq. (|30p . the linear 
conductance vanishes also in this case, G(V = 0) = 0. As 
we emphasized earlier, this result is due to an interference 
between the channels, i = ±. However, for a voltage in 
the range, < eV < T K , the Kondo effect in channel 
i = + is destroyed, while the Kondo effect is still well- 
developed in channel i = — . As a result, the destructive 
interference is suppressed as one turns on the voltage, 
and the conductance becomes non-zero. For even larger 
voltages, < eV, the Kondo effect is also destroyed in 
channel i = — , and the differential conductance starts to 
decrease with increasing voltage. Thus the result of the 
consecutive destruction of the two Kondo effects is the 
appearance of a peak at eV « T^ i 43 i 44 ' 56 This signature 
of the two-stage Kondo effect is clearly visible in the dif- 
ferential conductance, G(V), shown in Fig.Q3J There not 
only the central rcsonance-antiresonance structure can be 
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FIG. 12: Spectral functions for J/U = 0.2, T+/U = 0.4 
and Y-/U = 0.8, for different values of level splitting, A/J7. 
The upper panel shows the level-projected spectral functions, 
q±(uj), while in the lower panel the full spectral functions, 
qt(ui), are presented. Notice the breaking of electron-hole 
symmetry for any A / 0. 




0.6 0.8 
A/U 



FIG. 13: Linear conductance for J/U = 0.2, T+/U = 0.4 and 
Y-/U = 0.8, as obtained from numerical derivation of the 
current (ND), and from the Fermi liquid relations (FL). The 
inset shows the phase shifts 8±/tt extracted from the Green's 
functions at u — 0, and the occupation numbers, n±. 



seen, but also, a side-resonance at w w 2 J, which can be 
associated with the singlet excitations. 

As shown in Fig. [T51 for the lateral dot, increasing A 
gradually fills up the central dip, until the two resonances 
merge into a single ST-Kondo peak at around A w 2J, 
and finally split to two peaks located at eV « ±(A — 2 J). 
All these features are in excellent agreement with the 
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experimental results of van der Wiel et al.^ and Granger 
et al., (see Fig. 3e in Ref.|H). 

In Fig. [TS] we also present our results for the differen- 
tial conductance through a vertical dot. Contrary to the 
differential conductance of the lateral dot, these G(V) 
curves do not exhibit any remarkable difference with re- 
spect to the curves in Fig. [TU1 obtained for T + = r_. 
The main reason for this is that the perturbative ap- 
proach does not estimate correctly the Kondo tempera- 
ture: while for a non-perturbative calculation one would 
expect T£ <C for the couplings used, perturbation 
theory gives Kondo temperatures, Tj£, of the same order 
of magnitude. While this is sufficient to produce a visible 
effect in the differential conductance of a lateral dot, it 
is not sufficient to see the expected two-peak structure 
in the differential conductance of a vertical dot. In fact, 
for T£ <C Tg, the conductance should display two reso- 
nances on the top of each-other, a narrow one of width 
AV ~ T£ /e and a broad one of width AV~ ~ T^/e, both 
of height AG re 2e 2 /h. However, since X*j£ ~ Tj^ in our 
calculation, these two peaks merge into a single resonance 
of height Ae 2 /h on the triplet side of the transition. 

Our theoretical results for the vertical dots can be com- 
pared to the experimental curves of Sasaki et al£Q The 
evolution of the differential conductance on the singlet 
side of the transition and at the crossover point shows 
similarity with the experiments, but the data of Sasaki 
et al. do not display the three-peak structure expected 
on the triplet side. The reason for this is that, the experi- 
ments were performed at a finite temperature: the Kondo 
temperatures T K on the triplet side sink very quickly be- 
low the temperature as one gets farther away from the 
S-T degeneracy point, and therefore the central Kondo 
resonance was unobservable. 
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FIG. 14: Differential conductance of a lateral quantum dot for 
J/U = 0.2, V+/U = 0.4, and T-/U = 0.8, for level splitting 
A/U — 0. The side-peak can be attributed to the singlet ex- 
citation energy, 2 J, while the central resonance-antiresonance 
structure is associated with the two consecutive Kondo effects 
on the triplet side. 
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FIG. 15: Differential conductance, G(V), for J/U = 0.2, 
T+/U = 0.4 and T-/U = 0.8 for different level splittings, 
A/U for a lateral (a) and a vertical (b) quantum dot. The two 
consecutive Kondo resonances on the triplet side are clearly 
visible in the lateral dot (middle resonance-antiresonance 
structure), while the two Kondo resonances merge into a sin- 
gle peak in a vertical dot. 



IV. CONCLUSIONS 

In this paper, we studied the non-equilibrium singlet 
triplet transition in lateral and in vertical quantum dots. 
We described the quantum dots in terms of a simple two- 
level Anderson model, which included the Hund's rule 
coupling (J) as well as the charging energy of the dot 
({/), and we used non-equilibrium iterative perturbation 
theory to describe the transport through the dot. 

As a first step, we explored the range of validity of per- 
turbation theory and found that it breaks down at inter- 
mediate coupling strengths, where multiple solutions and 
hysteresis appear as an indication of spontaneous orbital 
polarization formation. Then we computed the equilib- 
rium spectral functions, determined the non-equilibrium 
Green's functions, and finally used the Meir-Wingreen 
formula to compute the differential conductance through 
the dot. 

To our great surprise, within its range of validity, 
this simple approach is able to describe essentially all 
experimentally-observed properties of the transition. In 
particular, for a lateral dot with level- dependent de- 
cay rates, T±, the differential conductance, G(V), dis- 
plays the correct resonance-antiresonance structure cor- 
responding to the two-stage Kondo effect on the triplet 
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side, while in the vicinity of the S-T transition these two 
resonances merge with the singlet excitation into a single 
Koiido resonance. On the singlet side of the transition we 
find a split Kondo resonance, and we recover the correct 
location of the side peaks, eV « ±(A — 2J). 

For a lateral dot, the linear conductance, Gu n (A) ex- 
hibits a broad maximum of height G max = 2e 2 /h around 
the transition, which is located at the correct level split- 
ting, A « 2J. On the other hand, for a vertical dot we 
find a maximal linear conductance of 4e 2 //i on the triplet 
side of the transition, which gradually crosses over to a 
small conductance on the singlet side. We also verified 
that our numerical calculations satisfy the Friedel sum 
rule, and that the linear conductance as computed by 
numerical differentiation of the Meir-Wingreen formula 
is identical to the one computed from the phase shifts 
using the Landaucr-Buttikcr formula (assuming a Fermi 
liquid ground state). 

The differential conductance of a vertical dot also be- 
haves in the correct way: A Kondo resonance (double res- 
onance) is found on the triplet side, which is suppressed 
and broadened around the transition, and finally splits 
into two resonances on the singlet side. 

We thus arrive at the remarkable and surprising con- 
clusion that the perturbative approach is able to describe 
this rather complex behavior, and is in rather good agree- 
ment with the experiments. It has however, some impor- 
tant limitations: (1) Its range of validity is limited to 
small and intermediate couplings, and (2) it is unable to 
account for the correct width of the Kondo resonance. 
Furthermore, though we checked that our solutions sat- 
isfy current conservation, simple perturbation theory is 
not a conserving approximation, and in general a lot of 
care must be taken to satisfy it. However, as we show 
in a consecutive paper^i the problem of the width of 
the Kondo resonance as well as that of current conser- 
vation can be overcome at the expense of getting some- 
what poorer resolution of the Hubbard peaks by using 
the method of fluctuation exchange approximation! 47 ' 60 
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Appendix A: Hybridized non-interacting Green's 
functions 

For completeness, here we enumerate the non- 
interacting Green's functions, g. They can be easily con- 
structed based on the Dyson equation 

g^H-g^H-SrH, (Al) 

where we used the matrix notation of the main text. The 
Green's functions, go, in Eq. (|A1[) denote the Green's 
functions of the isolated dot, with energies replaced by 
the effective energies, £i a , and its matrix elements are 
given by 

9oi„,vl' = s « £«'<W<W (w -£(„), (A2) 

with s K the Keldysh sign defined in the main text. The 
self-energy, Sp, denotes the self-energy coming from the 
hybridization of the dot and the leads, and its matrix ele- 
ments can be computed in terms of the tunneling matrix 
elements, U a , 

a=L,R 

^Ti<r,i'a' = — * 2tt &oo' ^2 Uati'a f a (co) , 
a=L,R 

^Tia,i'a' = -« 27T (W ^ t ia ti' a {f a {u) ~ 1) , 
a=L,R 

a=L,R 

with f a (ui) denoting the chemical potential-shifted Fermi 
functions of the leads. In the numerical calculations, we 
inverted numerically Eq. (|A1|) to obtain g. 
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